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We search for CP violation in a sample of 20000 Cabibbo-suppressed decays, D + — > K + KgTr + iv~ , 
and 30000 Cabibbo-favored decays, D+ -> K+K° s ty + tt . We use 520 fb 1 of data recorded by 
the BABAR detector at the PEP-II asymmetric-energy e + e~ collider operating at center of mass 
energies near 10.6 GeV. We search for CP violation in the difference between the T-odd asymmetries 
obtained using triple product correlations of the D + (Df) and D~ {Dj) decays, respectively. 
The T violation parameter values obtained are At{D + ) = ( — 12.0 ± 10.0 s ^ a ^ ± 4.6 g y S ^) x 10~ 3 

and Ar{Df) — ( — 13.6 ± 7.7 s t a t =t S-^syst) x 10 -3 , which are consistent with the standard model 
expectations. 

PACS numbers: 13.25.Ft, 11.30.Er 



In the standard model (SM) of particle physics, the vi- 
olation of the charge-conjugation and parity symmetries 
(CP) is introduced by the Kobayashi-Maskawa (KM) 
phase in the Cabibbo-Kobayashi-Maskawa quark mixing 
matrix [l| . The KM ansatz has been tested at high preci- 
sion in K and B decays, where the KM phase contributes 
to the quark transition amplitude at tree level. How- 
ever, further experimental efforts are needed in D meson 
decays, where CP-violating amplitudes are predicted to 
contribute to the observables at the 10 -3 level Q. 

The sensitivity to CP violation in D meson decays 
reached by the B factories is of the order of 5 x 10 -3 [3|- 



6]. Although this does not represent a measurement of 
SM CP violation, it provides a constraint on possible ef- 
fects beyond the SM. New physics models introduce CP 
violation in D meson decays through tree and one-loop 
diagrams. While predictions for CP violation in tree dia- 
grams are not different from those in the SM [0(1O~ 3 )], 
new physics in loop diagrams may enhance CP violation 
effects at the order of 10~ 2 Q- 

We report herein a search for CP violation in the de- 
cays D+ K + K Q s ir + %- and D+ ->■ K+K^tt+tt- us- 
ing T-odd correlations ||. We define a kinematic triple 
product that is odd under time reversal using the vector 
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momenta of the final state particles in the Dfs rest frame 



as 



Cn 



PK+ ■ (Ptt+ XPtt") 



(1) 



Under the assumption of CPT invariance, time-reversal 
(T) violation is equivalent to CP violation. 

We study the T-odd correlations by measuring the ob- 
servable expressed in Eq. (JT]) and then evaluating the 
asymmetry 



= r(c T > o) - r(c T < o) 
T ~ r(c T > o) + r(c T < o) ' 



(2) 



where V is the decay rate for the process under study. 
The observable defined in Eq. ([2j can have a nonzero 
value due to final state interactions, even if the weak 
phases are zero The T-odd asymmetry measured in 
the CP-conjugate decay process, At, is defined as 



, = r(-c T > o) - r(-c T < o) 
T = t{-c t > o) + r(-c T < o)' 



(3) 



where Ct = Vk- ' (p*-- x Pir+)- We can then construct 



At = -= {A T - A T ) , 



(4) 



which is an asymmetry that characterizes T violation in 
the weak decay process 1(J- 12 • 



At least four different particles are required in the fi- 
nal state so that the triple product may be defined using 
momentum vectors only [131 ] ■ The D meson decays suit- 
able for this analysis method are D + — > K + Kgir + ir~ , 
D+ -> K+K° s ir+ir-, and D° -> K +K~ ■*+■*-. The 
search for CP violation using T-odd correlations in 
D° — > K + K~it + it~ has recently been carried out by the 
BABAR Collaboration, and no evidence of CP violation 
has been observed 

Followin g th e suggestion by Bigi Q, the FOCUS Col- 
laboration [15j first applied this technique to a sample 
of approximately 500 reconstructed D + and D+ events, 
respectively. No evidence of CP violation was found. In 
the present analysis, we perform a similar measurement 
using approximately 2.1 x 10 4 D + and 3.0 x 10 4 Df me- 
son decay candidates. 

The analysis is based on a 520 fb _1 data sample 
recorded mostly at the T(45) peak and at center of 
mass (CM) energy 40 MeV below the resonance by the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
collider. Contributions to the data sample have been 
recorded near the T(3S) resonance (« 31 fb _1 ), and near 
the T(2S) resonance (s» 15 fb -1 ). In addition, two large 
samples of Monte Carlo (MC) simulated events have been 
analyzed. In these samples, the e + e~ — > cc produc- 
tion process is generated using Jetset7.4 1 161 . and the 
detector response is simulated by GEANT4 [17[. About 



1.1 x 10 9 generic e + e~ — > cc MC events, correspond- 
ing to 846 fb _1 , were generated to include the previ- 
ously measured intermediate resonances in the Df, de- 



W 

Df.X MC signal events 

(s) o 



cays, while 4.0 x 10 6 e + e 

(« 1025 fb _1 ), where X represents any system of charged 
and neutral particles compatible with the relevant con- 
servation laws, were generated in which the Df-. signal 
decays to K + Kgir + ir~ uniformly over the phase space. 
Both MC samples were processed using the same recon- 
struction and analysis chain as that used for real events. 

The BABAR. detector is described in detail else- 
where fl8j . We mention here only the subsystems used 
in the present analysis. Charged-particle tracks are de- 
tected, and their momenta measured, with a combination 
of a cylindrical drift chamber (DCH) and a silicon ver- 
tex tracker (SVT), both operating within the 1.5-T mag- 
netic field of a superconducting solenoid. The informa- 
tion from a ring-imaging Cherenkov detector, combined 
with specific energy-loss measurements in the SVT and 
DCH, provides identification of charged kaon and pion 
candidates. 

The D + and Df meson decay candidates are recon- 
structed in the production and decay sequence: 



e e 



XDts-.Dt, 

(s)> (s) 



K + K s ir + 7r-;K° s ^ir + TT-, (5) 



using the events with at least five charged particles. We 
7r + 7r~ candidates using a vertex and 

and re- 



reconstruct Kg 



kinematic fit with the Kg mass constraint \li 
quiring a x 2 probability greater than 0.1%. We accept 
only Kg candidates that decay at least 0.5 cm from the 
e + e~ interaction region (IR) and have a mass before the 
fit within 15MeV/c 2 of the nominal K s mass. The K s 
candidate is then combined with three charged-particle 
tracks with total net charge +1, to form a Df\ candi- 
date. We require the tracks to originate from a common 
vertex, and the \ 2 fit probability (Pi ) to be greater than 
0.1%. In order to improve discrimination between sig- 
nal and background, an additional fit is performed that 
constrains the three charged tracks to the IR. The x 2 
probability (Pa) of this fit is large for most of the back- 
ground events, whose tracks originate from the IR, while 
it is smaller for signal events, whose tracks originate 
from a secondary vertex detached from the IR, due to the 
measurable flight distance. Particle identification is 
applied to the three charged-particle tracks, and the pres- 
ence of a K + is required. Charged kaon identification has 
an average efficiency of 90% within the acceptance of the 
detector, and an average pion-to-kaon misidentification 
probability of 1.5%. We require the CM momentum of 
the Df^ candidate, p* , to be greater than 2.5 GeV/c. This 
requirement reduces the large combinatorial background 
from B decays, and improves the signal-to-background 
ratio significantly despite some loss in signal efficiency. 

We first study backgrounds from charm meson decay 
processes which yield the same event topology. 
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The decay D* + — >• 7r + £> produces a significant D 
peak in the KgK + ir~ mass distribution. A fit with 
a Gaussian signal function yields a mass resolution of 
<7 D a ^K° s K+-n- = 4-6 MeV/c 2 . Selecting Z?° candidates 
within zk'3(T D a^ K o K + n - of the D° mass, we observe a 
clear D* + peak in the distribution of the mass differ- 
ence Am — m(K + Kgir + ir~) — m(KgK + ir~). This con- 
tribution is reduced to a negligible level by requiring 
Am > 0.1465 GeV/c 2 . 

We also observe background from the decay D + — > 
-K'+.K'ij-K'tj, with one of the K s decaying into the bach- 
elor pions of Eq. ([5]). This contribution is removed by 
requiring the 7r + 7r~ invariant mass to lie outside a ±8.7 
MeV/ c 2 mass window around the nominal K s mass . 
We look for backgrounds from D + — > K s ir + Tr + ir~ de- 
cays by assigning a pion mass hypothesis to the kaon 
candidate. We observe a D + signal over a large back- 
ground. Simulation shows that this background pro- 
duces a broad structure in the high-mass region of the 
£>+ mass distribution. We also looked for background 
from A+ — > pK s TT + TT~ decay by assigning the proton 
mass to the kaon candidate. We see a signal over a 
large background. We find it impossible to remove the 
D + — > KgTr + ir + ir~ and A+ — > pKgTr + ir~ events with- 
out biasing our mass distributions. Our MC simulations, 
however, show that the presence of these backgrounds 
does not bias the extraction of the meson yields. 
As a further check, we select a high purity data sample 
(87.5 %) of D+ -> K° s it + tt + tt- decays and assign the K+ 
mass alternatively to both ir + . We compute the asym- 
metries on the resulting integrated distributions and find 
that they are all consistent with zero. A similar result is 
obtained when we perform the test on MC events. 

We divide the K + Kgir + ir~ mass spectrum into two 
regions in order to extract separately the D + and 
-D+ signal yields. For the former we require 1.81 < 
m(K + K° s iT + iT-) < 1.92GeV/c 2 , while for the latter we 
require 1.91 < m{K+ K°it + tt-) < 2.02GeV/c 2 . 

For further signal-to-background optimization, we ex- 
plore three variables: the CM momentum, p* , the differ- 
ence in probability, P\ — p2, and the signed transverse 
decay length, Lt = where d is the distance vector 

between the IR and the Dfs decay vertex in the trans- 
verse plane, and fir is the -D^ transverse momentum 
vector. Signal events are expected to be characterized 
by larger values of p* , due to the jetlike topology of 
e + e~ — > cc events, and larger values of Lt and P\ — P2, 
due to the measurable decay length. 

Figure[T]shows thep*, P1 — P2, and Lt distributions for 
signal and background in the D + and Df mass regions, 
respectively. The signal distributions are obtained from 
D + -> AT^7r + 7r + 7r~ and Df — > KgK~Tr + Tr + decays in 
data after background subtraction. These decay modes 
are kinematically similar to the signal modes, but have 



higher signal yields and better signal-to-background ra- 
tios. The background distributions in Fig.[T]are obtained 
from Dt-. — > K + KgTT + n~ sidebands in the mass distri- 
butions for data. 

The normalized probability distribution functions (V) 
of the three variables for signal and background are com- 
bined in a likelihood-ratio test 



(p*,P 1 -P 2 ,L T ) (6) 



to optimize the signal yields separately for D + and D+. 
The optimization of the cut is performed by maximizing 
the value of S/^/S + B, where S is the number of signal 
events and S + B is the total number of events in the 
signal region. The purity S/(S + B) of the peak improves 
from 11.2% to 51.4% and from 16.6% to 60.6% for D+ 
and Df , respectively. 

Figure shows the K + Kg7r + n~ mass spectra in the 
D + and Df regions before [(a) and (c)] and after [(b) 
and (d)] the likelihood-ratio test. For each region, the 
signal is described by the superposition of two Gaussian 
functions with a common mean value. The background 
is parametrized by a first-order polynomial in the D + 
region, and by a second-order polynomial in the Df re- 
gion. The fitted functions are superimposed on the data 
in Fig. and the fit residuals, shown above each distri- 
bution, are represented by Pull = (A data - A fit )/VA data . 
From these binned extended maximum likelihood fit, we 
extract the integrated yields N(D+ ) = 21210 ± 392 and 
N(D+) = 29791 ± 337 from the fits, where the uncer- 
tainties are statistical only. The mean value and width 
of the main Gaussian are [i D + = 1869.8 ± 0.1 MeV/ c 2 , 
a D+ = 3.76 ± 0.08 MeV/c 2 for D+ , and fi D + = 1969.0 ± 



0.1MeV/c 2 



3.67 ± 0.18 MeV/c 2 for Dj 



We next divide the data sample into four subsamples 
depending on charge and whether Ct (Ct) is greater 
or less than zero. We define 



N(D+ y C T >0) 
N(D+ y C T <0) 
N(D- sy C T >0) 
N(D- s)1 C T <0) 



2 
2 

N(D 



(*)' 



2 

N(D 



(l + 4r), 



(1 - A T ) , 



[l-An 



(7) 



and fit the corresponding mass spectra simultaneously to 
extract the yields and the values of the asymmetry pa- 
rameters At and At- In this fit, the shape parameters 
are shared among the four samples and are fitted together 
with the yields, N(D^) and N(D^ s - ) ), and the asymme- 
tries, At and At- The T- violating parameter At is then 
computed using Eq. 
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FIG. 1: Distributions of p* , Pi — Pi, and L T for D + (top panels) and Df (bottom panels) candidates. The distributions 
for signal and background are shown as solid and dot-dashed histograms, respectively. All distributions are normalized to 1. 
Signal distributions are extracted from D + — > Kg-K + ir + iT~ and D+ — ► KgK~ir + ir + for D + and Df decays, respectively. The 
background distributions are extracted from the — ► K + K%n + sidebands. 
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FIG. 2: The K + K'g-K^ ~tt~ mass spectrum in the D + mass 
region (a) before and (b) after the cut on likelihood ratio. 
Similar plots (c) and (d) are drawn for Df . The curves in 
(b) and (d) result from the fits described in the text. The 
distributions of the pull values are also shown. The x 2 / n dof 
values from the fits are 0.87 (£>+) and 0.95 (£>+). 



We validate the method by using the generic MC sam- 
ple. We find that the fit results for At, At, and the 
computed value of At are in good agreement with those 
in the simulation, both for D + and Df. 

All event selection criteria are determined before the 
final fit in order to avoid any potential bias. The true 
central values of At and At are masked by adding un- 
known random offsets. 

After removing the offsets, we measure the following 
asymmetries: 

A T (D+) = (+11.2 ± 14.1 stat ± 5.7 syst ) x lCT 3 , 



At(D-) = (+35.1 ± 14.3 stat ± 7.2 syst ) x 1CT 3 , (8) 



and 



A T (Dt) = (-99.2 ± 10.7 stat ± 8.3 syst ) x 1(T 3 , 
At(D-) = (-72.1 ± 10.9 stat ± 10.7 syst ) x lO" 3 . 



(9) 



We observe values of At and At which differ signifi- 
cantly from zero only for decay. This may indicate 
the presence of final-state-interaction effects for this de- 
cay process, perhaps as a result of the slightly different 
resonant substructure between D + and decay. For 
example, the K *°K* + final state can contribute only to 
D+ through a doubly Cabibbo-suppressed decay process. 
In the case of D + decay we find At and At to be con- 
sistent with zero, in contrast with the results of a similar 
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FIG. 3: Fits to the four D + K+ K° s tt + it' data subsamples. 
The pull values are shown above each mass distribution. The 
X 2 /n d of values from the fit are 1.07 {D + , C T > 0), 1.10 (D+ , 
C T < 0), 1.19 [D~, C T > 0), and 0.95 (D~ , C T < 0). 



FIG. 4: Fits to the four D+ -> 7\ + 7f|7r + 7r" data subsamples. 
The pull values are shown above each mass distribution. The 
X 2 /n do f values from the fit are 1.05 (£>+, CV > 0), 1.03 (D+ , 
C T < 0), 1.15 (D~, C T > 0), and 1.02 C T < 0). 



analysis performed on the corresponding D° decay sam- 
ple [3|: 

A T (D a ) = (-68.5 ± 7.3 stat ± 5.8 syst ) x 10" 3 , 
A T (D°) = (-70.5 ± 7.3 stat ± 3.9 syst ) x 10- 3 . (10) 

The fit results for the four data samples are shown in 
Figs. [3] and LU Using Eq. ((H) we obtain the T violation 
parameter values: 

A T (D+) = (-12.0 ± 10.0 stat ± 4.6 syst ) x 10- 3 (11) 

and 

A T (D+) = (-13.6 ± 7.7 stat ± 3.4 syst ) x lO" 3 . (12) 

For comparison, the value obtained for D° decay was [3} 

A T (D°) = (+1.0 ± 5.1 stat ± 4.4 syst ) x 10- 3 . (13) 

The sources of systematic uncertainty considered in 
this analysis are listed in Table HI and were derived as 
follows: 

1. We checked for possible asymmetries resulting from 
the detector response using large statistics signal 
MC samples in which the ^jt\ decays uniformly 
over phase space. These events are then weighted 
according to the resonant structures observed in 
the data (the resonances that contribute most are 
p° 7T+7T-, A"* K+tt-, and K*- -> K° s ir-). 



Small variations with respect to the generated val- 
ues are included in the evaluation of the systematic 
uncertainties. Using the same samples, we stud- 
ied the effect of the forward-backward asymmetry 
caused by the interference between the electromag- 
netic current amplitude e + e~ — > 7* — > cc and the 
weak neutral current amplitude e + e~ — > Z° — > cc. 
This interference results in a Df^/DV-. production 
asymmetry that varies linearly with the cosine of 
the quark production angle 6* , with respect to the 
e~ direction. Since the BABAR detector is asym- 
metric, the final and yields are not equal. 
To include this asymmetry in the MC samples, we 
weighted them for the cos 0* dependence measured 
in a previous analysis This study showed that 
the forward-backward asymmetry does not affect 
our measurements. 

2. We modified the likelihood-ratio selection criteria, 
and considered the observed deviations from the 
central parameter values as sources of systematic 
uncertainty. 

3. In order to check for final state radiation effects, we 
modified the fitting model by allowing the second 
Gaussian which describes the signal to have a free 
mean value. The background description was also 
modified by using higher order polynomials. 

4. The particle identification algorithms used to iden- 
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TABLE I: Systematic uncertainty evaluation for At, At and At in units of 10 3 for D + — ► K + Kgir + ir and Df —¥ 
K + K° S 7v+n-. 



Effect 


At (£>+) 


A T (£>+) 


At (D~) 


At (Df) 


At (Df) 


At (D-) 


(1) Reconstruction 


2.1 


2.8 


1.3 


0.7 


1.0 


1.3 


(2) Likelihood ratio 


1.1 


3.4 


5.6 


2.5 


7.8 


8.2 


(3) Fit model 


1.3 


1.1 


1.5 


0.1 


0.8 


0.7 


(4) Particle identification 


3.7 


3.3 


4.1 


2.2 


2.5 


6.7 


Total 


4.6 


5.7 


7.2 


3.4 


8.3 


10.7 



tify kaons and pions were modified to more strin- 
gent or looser conditions in different combinations. 

In the evaluation of the systematic uncertainty for each 
category, we keep the largest deviation from the reference 
value, and assume that the uncertainty is symmetric. It 
should be noted that the systematic uncertainty on At 
is not evaluated as the sum in quadrature of the errors 
on At and At- Instead, it is evaluated directly from the 
deviation of At resulting from the fits. This is why the 
error from the likelihood ratio or from particle identifi- 
cation is much smaller for At than would be expected 
from the uncertainties on At and At- 

In conclusion, we have searched for CP violation using 
T-odd correlations in high statistics samples of Cabibbo- 
supprcsscd D + — »■ K + Kgir + ir~ and Cabibbo-favored 
Df — > K + Kgir + ir~ decays. We obtained T-violating 
asymmetries consistent with zero for both D + and Df 
decays with sensitivities of « 1.0 % and « 0.8 %, re- 
spectively. We found that possible final-state-interaction 
effects in the K + KgTT + TT~ final state are larger for Df 
decay than for D + decay. 
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